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Abstract 



We present two results characterizing minimizers of the Chan-Esedoglu L 1 TV func- 
tional F(u) = J \Vu\dx + X J \u — f\dx; u, f : W 1 — > R. If we restrict to u = xs and 
/ = xa, e R n , the L lr TV functional reduces to E(E) = Per(S) + A|S A Cl\. We 
show that there is a minimizer £ such that its boundary 9S lies between the union of all 
balls of radius j contained in and the corresponding union of ^-balls in £l c . We also 
show that if a ball of radius ^ + e is almost contained in £1, a slightly smaller concentric 
ball can be added to S to get another minimizer. Finally, we comment on recent results 
Allard has obtained on L 1 TV minimizers and how these relate to our results. 

1 Introduction 

Much of the work in image analysis reduces to extracting and processing information from 
images. Much of that information is, in turn, carried by shapes present in the images. The 
methods for extracting information from images range broadly over stochastic, wavelet, and 
variational or PDE based methods. In the past five to ten years, the variational and related 
PDE methods have drawn a great deal of attention. 

In this paper we study one of these variational methods from a shape processing per- 
spective. More specifically, we establish new results concerning the properties of exact min- 
imizers for the rather new Chan-Esedoglu L lr TV functional. While this functional (which 
we now abreviate as simply the L 1 TV functional) is applicable to scalar functions on M n , 
we study the functional specialized to binary functions, i.e. binary images or shapes. 

The minimization of the L X TV functional, 



yields denoised images u that are smoothed yet close, in an L sense, to the measures image 
f (sometimes called the input image or noisy measurement). As is well known from studies 
of the Rudin-Osher-Fatemi total variation functional [7], 






(2) 
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the total variation term J \X7u\dx reduces oscillations while permitting sharp edges, some- 
thing that previous methods could not do or did very poorly. 

The change of the data fidelity term (J \u — f\ 2 dx in ([2])) to the L 1 term in the L 1 TV 
functional has the effect of making that functional much more natural from a geometric 
point of view. 

The L*TV functional was studied very carefully in a paper by Chan and Esedoglu [3]. 
(The discrete analog of the L 1 TV functional had been previously studied by Alliney [2] 
and Nikolova j6].) Chan and Esedoglu [3] show that for binary input images, there are also 
binary minimizers. More precisely, given any minimizer to ([I]) with binary input, almost 
every super-levelset is the support of a binary minimizer of the same functional. For binary 
input Xdi the functional can therefore be written as 

S* = argmin£(£) = Per(S) + A|S A U\, (3) 

where S, £* and O are the supports of the binary functions under study. Allard has recently 
submitted a paper |1J in which he uses very intricate geometric measure theory techniques 
to prove precise regularity results for minimizers of a class of functionals which includes 
the L 1 TV functional. We comment a bit more on Allard's work in the final section of the 
paper. 

Our results for minimizers of the L 1 TV functional can be viewed results on the reg- 
ularization of noisy shapes. The first result gives us a characterization of minimizers for 
the case in which the noise expresses itself as perturbations of the boundary. The second 
result characterizes the L 1 TV regularization of a binary images with measurement noise. 
In discrete images this corresponds to pixels flipping from to 1 or 1 to as driven by the 
noise process. 

Now a brief outline of the paper. In the next section we present the results for the case 
of S, £* and O in R 2 . This is the case most relevant for typical images. In the Section [3] we 
prepare for the proof of these results by introducing, in some detail, the notion of measure 
theoretic boundary, exterior and interior. This permits us to avoid the intricacies of the 
notion of reduced boundary. Next we prove the results for sets in R 2 (Section [4J . This 
section is the longest and most involved. In Section [5] we state the theorems for the case 
n > 2 noting a few modifications that must be made. Since all the hard parts of the proof 
for n > 2 are contained in the n = 2 case, we do not present the proof details. We close 
(Section [6]) with a brief discussion of our results and their relation to one of Allard's results. 

In what follows we represent minimizers of the L 1 TV functional ^ by S, dropping the 
superscript * used above. 

2 Main Results (n = 2) 

The two main results of this paper can be stated informally as follows. Define R = 2/A. 
For any ei, €2 > 0, 

(1) any ball of radius R completely contained in Q is also contained in E, and 
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(2) if a ball of radius R + e\ is almost contained in Q, then a concentric ball of radius 
R — €2 is completely contained in £. 

More precisely we have, 

Theorem 1. Let £1 be a bounded, measurable subset o/R 2 . Let S be any solution of ^5p. 
Assume that a ball Br of radius R lies completely in £1: Br C ft, . Then Br U T, is also a 
minimizer. Moreover, if Br C £l c , then (Br U S c ) c is also a minimizer. 

and, 

Theorem 2. Given f € (R, ^-R) and e E (0,1 — we can choose 5 = 5(R,r,e) > 

\B f \n\<6^ B (1 _ e)R c S. (4) 

Remark 1. Obvious analogs of these theorems hold in W 1 with modifications commented 
on in Section^ 

Remark 2. Theorem^ and the lower semicontinuity of the L1TV functional implies that 
there is a minimizer T, such that (J{-E>2 (x) Cl!}cS and (J{£?2 (x) C £l c } C X c . 

Remark 3. These theorems are close to optimal since the minimizer for f2 = B2 for 
arbitrarily small 77 > has unique minimizer S = 0. 



3 Measure Theoretic Boundary 

To simplify our analysis of the energy EiT,) = Per(S) + A|f2 A S|, we introduce measure 
theoretic boundary, interior, and exterior. 

Define Per(£) = J ' \^XT,\dx. We say a set in M n is a set of finite perimeter if Per (£) < 00. 
The structure theorem for sets of finite perimeter tells us that Per(S) = H n (d*Yi), where 
e?*X is the reduced boundary of S. The reduced boundary is rather complicated to define 
and difficult to manipulate. Instead, we use another theorem which asserts d*T, C d*X 
and H'^id^ - <9*£) = to conclude that Per(E) = fP-^&E), where <9*£ denotes the 
measure theoretic boundary of E. (See [5] Theorem 2, Section 5.7 and Lemma 1, Section 
5.8 for more details.) We now define measure theoretic boundary, interior, and exterior. 

Definition 1. A point x £ W 1 is in d*A, the measure theoretic boundary of A if 

C(B(x,r) n A) 

hmsup n > (5) 



and 



C(B(x,r)nA c ) 
limsup n > 0. (6) 

j — »o r 
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A point x S R n is in A\, the measure theoretic interior of A if 

£{B{x,r)DA c ) n 
hmsup n = 0. (7) 

while x £ M n is in A°, the measure theoretic exterior of A if 

C(B(x,r)nA) 

limsup v v \' - = (8) 

Lemma 1. Let A be a subset ofW 1 with finite perimeter. Then 

1. Per(A) = H n - l (d*A) 

2. R n = (A°) U (d*A) U (Al) and the three sets are pairwise disjoint. 

Proof: (1) As stated above this follows from [5] Theorem 2, Section 5.7 and Lemma 1, 
Section 5.8. (2) This follows directly from Definition [T] ■ 

Lemma 2. Suppose A and B be subsets ofW 1 . Then: 

1. if x £ A\ or x e B\ then x € (AU B)\. 

2. d*(AuB) c d*AUd*B 

3. d*A U d*B = (d*A n B\) U (d*A n B°) U (d*B n Al) U (d*B n Al) U (d*A n d*B) 

4. (1-3) immediately imply that d*(A U B) C {d*A n B°) U (d*B n A%) U (d*A n d*B) 

5. (d*A n B°) U (d*B n A%) C4(AU B) 

6. d*A c = d*A. 

7. (A C )Z = A\. 

Proof: The lemma follows in a straightforward manner from the definitions of measure 
theoretic boundary, interior and exterior. ■ 

Corollary 1. // H n ~ l (d*An d*B) = then 

1. H n -\d4A U B)) = H^id^A n B°) + H n - x {d*B n A°) 

2. H^id^A n 5)) = fl*- 1 ^*^ n sj) + H n -\d*B n 4) 
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Proof: (1): Lemma [2j (4-5) imply that 

H n ~ l {d*Af\ B°) + H n -\d*BnA°) (9) 

< H n -\d4AuB)) (10) 

< H n - 1 (d m AnB°)+H n - 1 (d*BnAZ) + H n - 1 (d m And*B) (11) 



and the conclusion follows. (2): This follows from (1), Lemma |2j(6)-(7) and the fact that 
A(~)B = {A C UB C ) C . 



Remark 4. Since d*A = (d*A n B\) U (d*A n d*B) U (<9*A n B°), the assumption that 
H n [d*A n d*B) = means we can, for the sake of measurement, consider d*A = (d^A n 
Bl) U (d*AC\B°). 

Remark 5. Now suppose that H n ~ 1 (d*A) < oo and B r is the ball of radius r centered at 
x G W 1 (we suppress the x). Note that d*B r = dB r . By the coarea formula, the set of r's 
such that H n ~ 1 (d^B r n d*A) > is at most countable. We conclude that the r's for which 
H n ~ 1 (d*B r n d*A) = are dense and have full measure in M.. For the rest of this section 
we assume that we have chosen r such that H n ~ (d*A fl d*B r ) = 0. 

Theorem 3. Suppose B r C O. Define E(E) = J \Vxx\dx + A / |£ s - fa\dx = Per(E) + 
A|EAO|. Then 

AE = E(H U B r ) - £(E) = -H n - l (d*X n (B r %) + H n ~ l {d,.B r n E°) ~ A|#r \ E|. 

Proof: Since Per(EU5 r ) = H n ~ l {dj:r\{B r )l)^H n ~ 1 {d, f B r r\T l %) andPer(E) = fP-^&En 
(.B r )i) + H™- 1 ^ n (B r )«) we get Per(E U (S r )) - Per(E) = -H n ~ l {d*Y, n (fl r )j) + 
ff n_1 (9*S r n E°). Noting that B r cn implies |(E U B r )AO| - |EAO| = |£ r \ E| finishes 
the proof. ■ 

The example A's and B's in Figure [T] illustrate why the above care is necessary. 

4 The Comparisons: proofs of Theorems [l] and [2] 

Recall that E denotes a minimizer of [3j If AS = E(B r U E) — i£(E) < then B r U E must 
also be a minimizer. 

Proof of Theorem [lj Computing AE for B r C O we get (for all but countably many r): 
AE= -iJ 1 (d*EnB4) + ff 1 (d*- B r-nE°)-A|.B r .\E| (12) 

= -tf^s n s r *) - tf 1 ^*^ n e;) + tf 1 ^*^ n e°) + H l (d*B r n sj) 

-A|S r \S| - A|5 r nE| +A|B r nE| (13) 
- J H" 1 (5*(B r nE))+ J H' 1 (^(S r ))-A| J B r | + A| J B r nE| (14) 
(H 1 (d^(B r )) — X\B r \) + (A|_B r n E| — H 1 (d*(B r n E))) . (15) 
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Figure 2: An illustration useful for computing AE when B r n C = 0. 
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This is illustrated in Figure [2] below. 



Now we choose a sequence of radii r, < R converging to R, for which (15) holds. 



Defining pi by Trpj = \B n f] E|, p* by 2-irp* = H l (d*(B ri n E)) and remembering that 
i |, we have that 

AE =(27rn-|7rrf) + (|7rpf-27rp*) (16) 
= 2^(1 -^) + 27r Pi (^-f). (17) 

Note that the isoperimetric inequality gives ^ > 1 for all i, that ^ < 1 for all i and 



that (1 — ^) ^j-too 0. The right hand side of (17) converges therefore to zero. Using the 
fact that AE is lower semicontinuous for sequences in L , (which follows from the lower 
semicontinuity of the BV seminorm), we conclude that AE(Br) < 0. We conclude that 
E U Br is also a minimizer. 

Finally we note that £q(E) = Per(E) + A|E A fi| = E Qa (E°) = Per(E c ) + A|E C A Q?\. 
From this we deduce that E minimizes Eq 44> E c minimizes Eqc Therefore, Br C £l c 
implies (E c U -Br) c is also a minimizer. ■ 

Proof of Theorem [2} In the case that B r ntt c ^ 0, 

ae = -H 1 (d*znB r i) + H 1 (d*B r nz°))-\\B r \z\ + 2\\B r nn c nz c \ (18) 
< -H 1 (d*znB r i) + H 1 (d*B r nz°)-\\B r \z\+2\\B r nn c \ (19) 

= -H 1 (d*ZnB r i)+H 1 (d*B r nZ°)- \\B r \Z\+2\\B r \n\ (20) 
= (^ 1 (^(5 r ))-A|5 r |) + (A|£ r n£| -^(^(BpnE))) +2A| J B r \0| (21) 
= (Per(B r )-A|S r |) + (A|B r .nS|-Per( J B r nS))+2A|B 7 .\n|. (22) 

This is illustrated in Figure [3] below. 

Lemma 3. |£r \ E| < 65 

Proof: Since we assume E is a minimizer, AE > 0. We will perturb with balls of radius 



r < r. Then, |£? r \ f2| < \B? \ £l\ '■= 5. These assumptions together with (22) and the 
isoperimetric inequality (Per(f? r D E) > 2^/n\B r D E|z) imply: 



i 



< AE < 2nr - Xirr 2 + X\B r n E| — 2y/n\B r n E|t + 2X5. (23) 
= A|5 r nS| -2V7r| J B r nE|5 + (2A<5 + 2vrr - Avrr 2 ) (24) 



f(0 = A^ 2 - 2^i + (2X5 + 2irr - Avrr 2 ) , (f=|S r nS|i) (25) 

2 /4<5 27rr 2 \ 2 
-£ 2 - 2-v/Trf + ( — + 2tit — (recalling R = —) . (26) 



i? V R R J X 

In view of the last inequality, we describe values of £ for which /(£) > 0. For a given r, the 
zeros of /(£) are at: 



t±(r) = \!~ ± a/ 7 ^ 2 +(7rr(r-fl))-25. (27) 
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d,B r n s° 





3 = \B r n!! c ns c | 

\B r nU c \ = 3 + 4 

4 = \B r n n c n s| 



a*B r n El 



Figure 3: An illustration useful for computing AE when B r f] £l c ^ 



Thus, for all r < f, we have: 

Either \B r n E|5 < £_(r), or |J3 r n S|i > £ + (r). 
If we take r = f > R, then 2nr(r — R) > 0, and assuming 
Condition 1. 



implies £-( r ) < 0. This implies 



/7 

Since r = r < ^~R we get 



7rr 

5< — (r-i2) 



l-Bf nE| >£+(r) >vri? 2 . 



3vrii 2 vri? 2 



B fl nE > LB r n£ — > 



Now we consider £±(R): 
Assuming 



irR 2 irR 2 

U(R) = \I — ±\I—-2S- 



(28) 



(29) 



(30) 



(31) 



(32) 
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Condition 2. 

S<^, (33) 



we get that £±(R) are real and distinct. Since 



we conclude that 
Computing, we get 



irR 2 

£(R) < — < \B R nX\, (34) 



e+<\B R nE\. (35) 




ttR 2 ttR 2 



>- + — S<_^) (3S) 

= ttR 2 - (2 + 4a)S. (39) 



Choosing a = 1, and noting that Condition [2] then implies the assumption in 38 is satisfied, 
we get 

\B R nE\ >£+ > tt-R 2 - 65. (40) 

This gives 

|£ fl \E|<6<5 (41) 

as advertised. ■ 



Remark 6. What if either r or R are radii such that (22) (and therefore (21)) does not 
hold? We can simply choose another f < f arbitrarily close to f, for which \22^ does hold. 
The 5 = \Bf\Ct\ will be no greater than, and arbitrarily close to, 5. As we will see, the only 
conditions on 5 that are not functions of R and e are those in Condition [IJ Therefore, if 
we replace Condition^ with 

5<^(f-R) (42) 

we know that the delta chosen for any f will permit us to arrive at the conclusions of this 
lemma, even in cases where we have to perturb r. Next we choose a sequence of r{ > R 
converging monotonically to R for which the inequality does work. Equation \31\) is still 



valid if we replace B R with B ri . Equation (32) can be slightly modified using \21^ to 



^) = J^±J^-25 t . (43) 
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where the 5i < 5 and 5i — > 5 as i — > oo. Now, simply repeating the derivation in lines (36) 



to (39), gives 

\B R n E| = lim \B n n E| > ttR 2 -Q5= lim ttR 2 - 65i (44) 

i— »oo i— >oo 

Now we continue with the proof of Theorem [2j Computing (again and less optimally, 
but sufficiently for our purposes) the change in energy when we add a ball B r to S for 
r G (0, /?), we get 

AE = E(T,U B r ) - E(E) (45) 

< -fl' 1 (a,EnB r i) + £r 1 (a l ,B r ns;) + A|B r \s| (46) 

= -Per(S ;J B r ) + J ff 1 (a^ r nS°) + A|B r \S|. (47) 
By the coarea formula and properties of the measure theoretic exterior, 

\B r \Z\ = \B r nX°\= [ H^d^nE ,)^. (48) 

Jo 

By the relative isoperimetric inequality applied in the ball B t {xq), 

Per(S ;J B r ) > CnunjlBrXSISjlSnBrlH . (49) 

Assuming 

Condition 3. Q5 < \-kR 2 

implies |-Br \ S| < ^ttR 2 . Assuming r > ^= implies that |B r \ S| < |E D -B r | and 
consequently 

Per(£;£ r ) > C|S P \E|2. (50) 

This gives a condition on e: 
Condition 4. e < 1 7=. 

v2 



Define w(r) := |5 r \ E|. By differentiating (48) with respect to r, and using (50) we see 



that the inequality concerning the change in energy given in ( 45 ) can be written as 



E(H U B r ) - E(E) < Xv(r) - Cy/v^r) + v\r). (51) 

We will use the differential expression on the right to show that the change in energy on 
the left has to be negative for some r close to R. 

Remark 7. Note that by choosing 5 small enough, we can make v(r) arbitrarily small and 
obtain Xv(r) — Cy/v(r) < 0; if the right hand side is positive then we have v'(r) > 0. This 
in turn means that v{r) decreases as r gets smaller. We exploit this to force the right hand 
side to zero. 
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Lemma 1. v'(r) — C^v{r) + Xv(r) < for a set of ' r £ ((1 — e)R, R) with positive measure. 

Proof of lemma: Assume 

v'(r) - C^/vJr) + \v{r) > for a.e. r £ ((1 - e)R, R) , (52) 

otherwise we are done. Let w(s) := e~ Xs v(R — s). Then (52) turns into 

w'(s) + Ce^ \/w{s) < for a.e. s £ (0, eR) . (53) 

with the initial condition w(0) = \B R \ E| and w(s) > 0. Solutions of this differential 
inequality can be bounded from above by solutions of the following differential equality: 



— / — AS i—— 

w = —Ce^^yw. 



w(0) = \B R \ E| and w > 0. 

The solution is 



(54) 



y^Msj = max (V j (e~^ s - l) + y/\B R \ E\\ 
= max (o, ™ (e"S - l) + y/\B R \ EM 



Therefore if \B R \ E| < 65 and 

Condition 5. 65 < a, where a is any solution to 



CRf eR \ r- CRf f 

' e~^ -1 )+^a= — e" £ -l + y 7 " < (55) 



2 I' J v 2 V 

i.e., we have 

C 2 P 2 

5<^(l-0 2 , (56) 

then we have a set of r with positive measure in ((1 — e)R,R) such that v(r) = and 
v'(r) = 0. ■ 

This lemma immediately implies that for some r £ ((1 — e)R,R), B r U E is also a 
minimizer. ■ 



5 The Case n > 2 

The analogs for Theorems [T] and [2] in W 1 are: 

Theorem 4. Let Q be a bounded, measurable subset o/R n . Let E 6e any solution of 
Assume that a ball B R , R = ? o/ radius R lies completely in £1: B R C f2 . T/ien 1?^ U E is 
afeo a minimizer. Moreover, if B R C i/ien {B R U E c ) c is a/so a minimizer. 
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Theorem 5. Given 



and 



f G \R, 2 



n 



n 



R 



e G 0,1 



we can choose 5 > suc/i i/iai 



|J3^ \ n| < S B ( 



l-e)R 



U £ is a/so a minimizer. 



(57) 



We do not present the proofs, since they are very similar to the n = 2 case. In particular, 
making the replacement R = j — > R = j enables us to use the proof of Theorem [l] with 
obvious modifications, to obtain Theorem |4j Likewise, we can use the proof of Theorem [2] 
to prove Theorem [5] with modifications noted below. 



n 
X' 



(1) Again, R= £ -> R 

(2) We define ^ = \B r n E|«. Let ce n be the volume of the ball with unit radius in 



The polynomial in (26) then gets replaced by 



/r + {-R- + nanr r 

Since we are interested in the roots of this polynomial, we look at 



n-l 



na n r 



aZR-t) =25 + a n r n - 1 {R 



(3) We replace the right hand side of ( |51[ ) with 

n 1 

\v{r) — Cv(r)~"^ + v (r) 

which gives us 

— As n-l 

w'(s) + Ce~ (w(s))^ < 



(58) 
(59) 

(60) 
(61) 



in place of (53). 



6 Discussion 

As mentioned in the introduction, Allard [1] has recently produced an extensive study of 
the regularity of minimizers for a class of functionals including the L 1 TV functional. In this 
work he uses geometric measure theory techniques originally developed to address minimal 
surface problems. As a result, his n cannot exceed 7. In our work we have used simpler pieces 
of geometric measure theory, specifically the structure theory for sets of finite perimeter. 
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The weaker regularity results we use - simply what one gets from £ having finite perimeter 
- are not limited to n < 7. 

Allard proves that the total mean curvature of minimizers is bounded by A, as suggested 
by a naive calculation with the formal Euler-Lagrange equation. For spheres, this corre- 
sponds to a radius of curvature of ^r— • In our work we find that spheres (balls) of radii j 
play a critical role. Such a sphere has a total mean curvature of ^F— A. This second, bigger 
radius characterizes the global nature of the minimizers. Indeed if one can contain £1 in a 
ball of radius j — e, where < e, then the unique solution is the empty set. This follows 
from a monotonicity result proved by Yin in [9J. It also follows from monotonicity results 
in Allard's paper pQ. 

Another previous work that needs to be mentioned is the work of Italo Tamanini and 
collaborators (see [B] and references). Instead of using knowledge of f2 to deduce properties 
of the minimizer X, they use weaker properties of £ to establish stonger properties of the 
same S. (In particular, if a ball of a particular radius is almost contained in £ then the 
ball with half the radius and same center is completely contained in S. This is similar 
to our Thereom[2j) These types of regularity properties of minimizers are not very useful 
for computation or in the establishment of minimizer properties based only on realistically 
obtainable knowledge. 

There are numerous potential directions in which to advance to these results and the 
results reported in [U|9]. Generalization to anisotropic energies (see |3j for example), the 
construction of hybrid analytic- numerical algorithms for L 1 TV minimization, and the ex- 
ploitation and analysis of the scale decomposition properties of the L 1 TV functional are 
three that come easily to mind. 

7 Acknowledgements 

It gives me pleasure to acknowledge significant conversations with Selim Esedoglu who 
contributed substantial suggestions for improvment as well as suggesting the initial idea 
that started this paper. I also acknowledge the pleasant benefit of discussions with Bill 
Allard. The research was made possible by funding from Los Alamos National Laboratory 
and the Department of Energy. I would also like to acknowledge the hospitality of the 
Institute for Pure and Applied Analysis at UCLA where I was able to work uninterrupted 
on the research reported here. 

References 

[1] William K. Allard. Total variation regularization for image denoising; I. Geometric 
Theory. To appear in SIAM Journal on Mathematical Analysis, 2007. 

[2] S. Alliney. A property of the minimum vectors of a regularizing functional defined by 
means of the absolute norm. IEEE Trans. Signal Process., 45:913-917, 1997. 



IS 



Kevin R. Vixie 



Some Properties of CE-L TV Minimizers 



[3] Tony F. Chan and Selim Esedoglu. Aspects of total variation regularized L 1 function 
approximation. SI AM J. Appl. Math., 65(5):1817-1837, 2005. 

[4] Selim Esedoglu and Stanley J. Osher. Decomposition of images by the anisotropic 
Rudin-Osher-Fatemi model. Commun. Pure Appl. Math., 57:1609-1626, 2004. 

[5] Lawrence C. Evans and Ronald F. Gariepy. Measure Theory and Fine Properties of 
Functions. Studies in Advanced Mathematics. CRC Press, 1992. ISBN 0-8493-7157-0. 

[6] Mila Nikolova. Minimizers of cost-functions involving nonsmooth data-fidelity terms. 
SI AM J. Numer. Anal, 40:965-994, 2003. 

[7] Leonid Rudin, Stanley Osher, and Emad Fatemi. Nonlinear total variation based noise 
removal algorithms. Physica D, 60(1-4) :259-268, November 1992. 

[8] Italo Tamanini and Giuseppe Congedo. Optimal Segmentation of Unbounded Functions. 
Rend. Sem. Mat. Univ. Padova, 95:153-174, 1996. 

[9] Wotao Yin, Donald Goldfarb, and Stanley Osher. Image cartoon-texture decomposition 
and feature selection using the total variation regularized L 1 functional. Submitted to 
SIAM MMS, 2006. UCLA CAM Tech Report 05-47. 



